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Kliem and Torok Reply: The Comment by Chen 
[H on our investigation of the torus instabihty (TI) of 
arched magnetic flux ropes (hereafter KT 0) states that 
he has previously solved a more general equation than our 
Eq. (1) — Eq. (7) in Ref. [1| — (including the effects of non- 
vanishing plasma beta), using less restrictive assumptions 
(not assuming constant aspect ratio so that "the 

'torus instability' is a limiting case of Ref." In addi- 
tion, it is claimed that {i) the absence in KT of a term 
that introduces a spatial scale (the footpoint distance 
Si of an expanding rope whose ends are anchored in a 
rigid surface) makes our model conflict with the presence 
of such a scale in the observations of coronal mass ejec- 
tions (CMEs) and {ii) the assumption in KT's Eqs. (8-10) 
that the ring current / tends to be conserved as a con- 
sequence of the rope's footpoint anchoring is mathemat- 
ically inconsistent with the ideal-MHD assumption that 
the magnetic flux \E' enclosed by the rope is conserved. 

The TI is a property of the tokamak equilibrium Q 
in the low-beta limit, which is the relevant case in the 
coronal source regions of CMEs. For /? < 1, the current 
ring can attain equilibrium only in the presence of an 
external poloidal field Box tJ]- The instability occurs if 
-Box decreases sufficiently steeply with increasing major 
radius R, with a decay index n = — dlnScx/c^ln-R ^ 
3/2 0, [g. In Ref. the simplification Sox = was 
adopted in the starting equations of both, the linear and 
the nonlinear analysis (Eqs. [7] and [32], respectively), 
thus excluding the TI. (Moreover, since this requires f3 > 
1 for equilibrium 4], the results in Ref. Q are largely 
irrelevant for the explanation of CMEs.) 

Unlike Refs. [1, [E 0]' replaced the integration of 
the force equation in the b direction by the assumption 
of self-similar expansion, R/b = const. This is required in 
order to permit an analytical description (our Eqs. [4-10], 
which have no analogue in Refs. iiSi), and is also 
a very reasonable simplification, suggested by the obser- 
vations and justified by the facts that the hoop force de- 
pends only logarithmically on R/b and that our resulting 
description yields qualitative and quantitative agreement 
with essential CME properties ^^jjA, ilfl] . Chen adopted 
this assumption recently as well [8|, based on the expe- 
riences made in fitting his model to CME observations. 
In integrating the equation for b(t) in Refs. [1, @, the 
energy equation was replaced by the polytropic assump- 
tion, with the index 7 « 1.2 estimated from comparison 
with observations; this is a level of approximation com- 
parable to ours. Our assumption R/b = const obviously 
implies (Pb/dt^ oc (PR/dt"^, not d'^b/dt^ ~ as stated in 
the Comment. 

It is well known that the "standard" MHD equations 
(e.g., 0]) do not contain an intrinsic length scale. Di- 
mensional lengths enter by prescribing initial or bound- 
ary conditions at the application stage, or by specifying 
the treatment such that these conditions enter already 
explicitly at an intermediate stage. When our scale-free 
Eq. (4) is applied to describe CMEs, the footpoint an- 


choring of the flux rope, which was not explicitly included 
in KT, introduces the condition of nearly semicircular 
flux rope shape at TI onset (which is supported by ob- 
servations [nj and was suggested in Ref. [t^ as well); 
i.e., i?o ~ Si/2. With the peak acceleration occurring in 
the range R ^ (1.5-2)i?o ~ (0.75-l)S'f in the practically 
most relevant range of n included in Fig. 1 in KT, our 
description reproduces the observational result of Fig. 1 
in the Comment to a very reasonable approximation, the 
more so if the spread in the position of observed ejecta 
relative to the magnetic rope axis is taken into account, 
and it does so in a manner that is more general than 
the use of a problem-specific inductance in Refs. [1, 0, S| 
(which can easily be incorporated in our theory). Con- 
trary to a statement in the Comment, we did not suggest 
any connection between I{t) and the scale Sj. 

I and ^E* can be simultaneously constant in a certain 
range of R. This is obtained by using these two condi- 
tions joint with our general ansatz Bcx{R) = BR~" in 
Eq. (2) and elementary algebra to solve it for b/R, given 
in Eq. (10) and plotted in Fig. 3 of KT. Using constant 
/, our Eq. (8) follows exactly from our Eq. (1). The 
false conclusion in the Comment follows from the inap- 
propriate simplification L oc R, neglecting the ln(i?/6) 
dependence, which removes b from the equation. 

Our estimate of the instability threshold ricr for con- 
stant / (Eq. [9]) includes an inconsistency because R/b = 
const was used in addition to constant I and How- 
ever, (i) we have expressed that Eqs. (8-9) represent a 
limiting case, included to demonstrate the direction of 
the effect of constant /, (ii) this estimate does not play 
any further role in KT, and (Hi) the inconsistency of this 
approximation is made apparent immediately following 
Eq. (9) and Fig. 2 in KT. 


B. KHem^ and T. Torok^ 

^Kiepenheuer Institute, Freiburg, Germany 
^LESIA, Observatoire de Paris, CNRS, France 

Received: 23 October 2006, revised 2 Aug/21 Sep 2007 
PACS numbers: 52.35.Py 52.30.-q, 96.60.ph, 96.60.qf 


[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 


J. Chen, Phys. Rev. Lett. 99, 099501 (2007). 

B. Kliem, T. Torok, Phys. Rev. Lett. 96, 255002 (2006). 
J. Chen, Astrophys. J. 338, 453 (1989). 

V. D. Shafranov, Rev. Plasma Phys. 2, 103 (1966). 

G. Bateman, MHD Instabilities (MIT, Cambridge, 1978). 

J. Chen, J. Geophys. Res. 101, 27499 (1996). 

J. Chen and J. KraU, J. Geophys. Res. 108, 1410 (2003). 

J. Chen et al., Astrophys. J. 649, 452 (2006). 

T. Torok, B. Kliem, Astron. Nachr. 328, 743 (2007). 

C. J. Schrijver, C. Elmore, B. Kliem, T. Torok, A. M. 
Title, Astrophys. J. subm. (2007). 

B. Vrsnak et al., Solar Phys. 136, 151 (1991). 


